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Abstract 

A randomly pinned elastic medium in two dimensions is modeled 
by a disordered fully-packed loop model. The energetics of disorder- 
induced dislocations is studied using exact and polynomial algorithms 
from combinatorial optimization. Dislocations are found to become 
unbound at large scale, and the elastic phase is thus unstable giving 
evidence for the absence of a Bragg glass in two dimensions. 

Randomly pinned elastic media are used to model various condensed- 
matter systems with quenched disorder including the vortex phase of dirty 
type-II sup er conductor s|l]. Much analytical progress on these systems has 
been made within the elastic approximation where dislocations are excluded 
by fiat. The intriguing possibility of spontaneous formation of disorder- 
induced dislocations (pairs and loops respectively in two and three dimen- 
sions) at large scale, however, remains a challenging question ||. To address 
this issue at zero temperature requires a detailed understanding of the ener- 
getics of dislocations in terms of their elastic-energy cost and disorder-energy 
gain in the ground stateBG. 

In recent years, we have witnessed a fruitful exploration of novel al- 
gorithms in for complex disordered systems whose ground state itself is 
dominated by random disorder. One class of these efficient algorithms is 
based on network flow optimization. It includes the min-cost-flow, max-flow 
and matching algorithms which compute the exact ground state in time 
that grows only polynomialy in the system size, an attibution of great prac- 
tical importance. Some recent applications include studies of the rough- 
ness and topography of random manifolds and 2d random elastic media 
by max- flow ||], matching [|j and min-cost-flow algorithms ||, the sensitiv- 
ity exponents of the random-field-Ising model by max-flow algorithms |J, 
the domain- wall energy in the gauge glass by min-cost-flow algorithms B, 
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the droplet excitations in disordered systems by matching algorithms^, 
and the critical exponents of 2d generic rigidity percolation by matching 
algorithms |J. 

In this article we briefly review our recent work ||i"o| on dislocations in 2d 
randomly pinned elastic media and show how the energetics of dislocation 
pairs can be studied numerically by applying these polynomial algorithms to 
a 2d lattice model. The essential ingredients required of such a 2d discrete 
model would be: (1) its large-scale fluctuations are described by an elastic 
Hamiltonian with a quenched random potential that reflects the periodicity 
intrinsic to any elastic medium; (2) dislocations can be "conveniently" gen- 
erated; and (3) its ground-state energies with and without dislocations are 
amenable to exact numerical computations by these polynomial algorithms. 



Models Fortunately, recent works of Henley, Kondev and their co-workers 
provided us with a large class of just such models whose degrees of freedom 
are described in terms of colors, tilings (dimer) and loops[ll|, precisely the 
natural language for considering network flow optimizations. More impor- 
tantly, these models all permit a solid-on-solid (SOS) representation whose 
large-scale height fluctuations are governed by a few elastic constants and a 
locking potential that is periodic in heights. As an illustration, we consider 
here a fully-packed loop (FPL) model defined on a honeycomb lattice. All 
configurations of occupied bonds which form closed loops and cover every site 
exactly once are allowed, as in the example of Fig. 1(a). The corresponding 
SOS surface is a (lll)-interface of a simple cubic lattice constructed as fol- 
lows. Define integer heights at the centers of the hexagons of this honeycomb 
lattice then orient all bonds of the resulting triangular lattice connecting the 
centers such that elementary triangles pointing upward are circled clockwise; 
assign +1 to the difference of neighboring heights along the oriented bonds 
if a loop is crossed and —2 otherwise. This yields single-valued heights up 
to an overall constant. 

We introcude quenched disorder via random bond weights on the hon- 
eycomb lattice, chosen independently and uniformly from integers in the 
interval [—w, w] with w = 500. The total energy is the sum of the bond 
weights along all loops and strings. The FPL model is shown to be equivalent 
to an array of fluxlines confined in a planeQ with the heights corresponding 
to the displacement fields of the fluxlines. The SOS surface described above 
can be viewed as an elastic surface embedded in a 3d random potential that 
is periodic in heights modulo 3 since the smallest "step" of the surface is 
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Figure 1: The FPL model with periodic boundary conditions. The ground 
states with and without a dislocation pair for one realization of random 
bond weights are displayed in (b) and (a) respectively. The dislocations 
(solid dots) in (b) are connected by an open string (thick line) among the 
loops. The relevant physical object is, however, the domain wall which is 
induced by the dislocations as shown in (c). This domain wall represents 
the line of all bond differences between the ground states (a) and (b). 



three[ll]. The coarse-grained effective Hamiltonian becomes [|T0|, 11, 12\ 
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where the random bond weights enter as random phase 7(r). Note also both 
K and u depend on the disorder strength w since it is the only energy scale 
in the problem. This is the well studied model for charge-density waves 
(CDW)||. 

Dislocations are added to the FPL model by "violating" the fully-packed 
constraint. One dislocation pair is an open string in an otherwise fully- 
packed system as shown in Fig. 1(b). The height change along any path 
encircling one end of the string is the Burgers charge ±3 of a dislocation 
so that the heights become multi-valued. Note that the configurations with 
and without a dislocation pair only differ along a domain "wall" , as shown 
in Fig. 1(c). Dislocations with higher Burgers charges ±6 can also be created 
by introducing holes instead of strings. 



Algorithms It turns out that the ground states of the disordered FPL 
model can be obtained via polynomial algorithms. A general description of 
such an optimization problem is given by the so-called linear programming 
which is to identify a set of variables minimizing a linear objective function 
subject to a set of linear constraints. Most physical problems are restricted 
to integer-valued variables. Such an integer optimization problem in gen- 
eral is nondeterministic polynomial (NP) which implies that polynomial and 
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exact algorithms are unlikely to be found. However, for a special class of 
problems where the linear constraints other than the upper and lower bounds 
on the variables can be interpreted as "flow conservation" at the nodes of 
a graph while the variables are identified with the flows on the edges of the 
graph, the optimization problem is polynomial. Recent applications men- 
tioned above fall into this class. Since most textbooks [[lj] contain details 
on the proof of this result and existing C++ codes for these polynomial 
algorithms can be found in the LEDA library[15], we shall only discuss how 
to transform the search for the ground states into an integer min-cost-flow 
problem on a suitably designed graph. The min-cost-flow problem is to find 
the flow pattern of minimum total cost for sending a specific amount of flow 
from a given node s to another given node t in a graph G in which the flow 
x on every edge has an upper bound Ub and a lower bound lb (h < x < u\>) 
as well as a unit cost c. The total cost is of course given by summing cx 
over all edges in G. 

Suppose that the bipartite honeycomb lattice with periodic boundary 
conditions contains 2N sites which we divide into two sublattices of N A- 
sites and N B-sites. We can construct a graph G as follows. In addition 
to all sites and bonds of the honeycomb lattice, this graph contains two 
extra sites, denoted as s (the source) and t (the sink), and extra 2N bonds 
(the leads). All bonds of the honeycomb lattice are directed from A-sites 
to B-sites with lb = 0, Ub = 1, and c being the corresponding random bond 
weight, while the remaining N in-leads are directed from s to A-sites and 
N out-leads from B-sites to t with l u = 1, ii& = 2, and c = for all 2N 
leads. Therefore, the ground state energy of a loop configuration with or 
without defects is equivalent to the minimum-cost flow if loops and strings 
are identified with bonds on the honeycomb lattice that have flow (note 
that the flow value on these bonds must be either zero or unity). Simple 
inspection shows that this identification can indeed be made with the above 
choice of bounds if the amount of flow sustained between s and t, with flow 
conservation on all other nodes, is between N and 2N units. 

Given the amount of flow sustained, the minimum-cost-flow algorithm 
establishes the flow pattern of the minimum cost. Various interesting phys- 
ical situations can be simulated by simple variations. 2N units of flow, for 
example, lead to the ground state of fully-packed loops (no dislocations). 
2N — 1 units of flow, on the other hand, give the ground state with one 
dislocation pair without a priori fixing the pair location. Keeping 2N — 1 
units of flow while changing Ub of a particular in-lead and out-lead from 2 
to 1 simulates a fixed pair of dislocations with the Burgers charges ±3. If 
using 2A — 2 units of flow instead and changing Ub of a particular in-lead 
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Figure 2: Energetics of a dislocation pair. Diamond and square symbols 
denote the defect energy for a pair of fixed dislocations with the Burgers 
charges ±3 and ±6 respectively. Solid lines are linear fits. Data denoted by 
circles are the defect energy E^ for a pair of optimized dislocations with the 
Burgers charges ±3. 

and out-lead from 2 to 0, we obtain a pair of dislocations (holes) with the 
Burgers charges ±6 at fixed locations. Clearly, dislocations of any desired 
density can be achieved by suitably varying the flow between N and 2N 
units. Moreover, introducing another extra link from t back to s with a 
negative unit cost — E c allows us to determine the optimal amount of flow 
sustained (thus the optimal dislocation density) with E c being the core en- 
ergy. This last simple variation results in the min-cost-circulation problem 
in network flow optimization. 

Numerical results For a given disorder realization, two ground-state en- 
ergies, E\ and Eq, were obtained respectively for cases with and without 
dislocations. The defect energy Ed = E\ — Eq was then determined. Various 
L x L sample sizes with L = 12, 24, 48, 96, 192, and 384 (480 for optimized 
defects) were simulated with at least 10 disorder averages for each size. 

We first describe our results for a single dislocation pair where the core 
energy E c is set to zero. The elastic constant of an elastic medium can be 
measured in various ways by observing its response to perturbations. Here 
we perturb the system with a fixed dislocation pair (large-scale topological 
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Figure 3: Optimal dislocation density p as a function of the core energy 
E c , Shown in the inset are the elastic-energy costs (E^ = E e i a ) to a pair of 
fixed dislocations with E c = and the Burgers charges ±3 injected into a 
state with already optimal number of pre-existing dislocations (square) as 
well as into a state with no pre-existing dislocations (circle). Data denoted 
by circles are the same as those denoted by diamonds in Fig. 2 and are 
shown here for comparison. 

excitation). The defect energy Ed in this case is the elastic energy cost E e i a 
which according to the elastic theory should scale as E e [ a ~ Kb 2 /2i: ln(L). 
This is indeed consistent with our numerical results shown in Fig. 2, and the 
elastic constant K is found to be 126(2) and 125(1) from dislocations pairs 
with the the Burgers charges ±3 and ±6 respectively. When the dislocation 
pair is allowed to be placed optimally, Ed also contains the disorder energy 
gain Edis in addition to E e i a , i.e., Ed = E e i a + Edis- As shown clearly in 
Fig. 2, Edis dominates over E e \ a resulting in the negative Ed, and moreover, 
Edis drops faster than ln(L). Detailed analysis showed that the numerical 
results are consistent with the theoretical prediction E^ a ~ — In 3 / 2 (X) ||10|| , a 
result independent of the disorder strength w. Therefore the elastic phase of 
large systems is unstable to dislocation pairs. With no restrictions on their 
number, dislocations will proliferate thereby driving the elastic constant K 
to zero. 

We now discuss our results on multiple dislocations which are summa- 
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rized in Fig. 3. It is indeed clear from the inset that the elastic energy 
cost E e i a for introducing a pair of fixed dislocations into the state where the 
number of dislocations is already optimal is independent of the separation of 
the fixed pair implying a zero elastic constant K. This is consistent with the 
result on a related model) 16]. The relation between the optimal dislocation 



density p and the core energy E c is, however, found to be 

with a = 0.74(3). This exponent remains elusive to us at the present. 



Conclusion and Outlook In conclusion, we studied the energetics of 
dislocation pairs in a 2d random elastic medium by applying polynomial 
algorithms to 2d disordered FPL model and found the elastic phase is un- 
stable against the proliferation of dislocations, and thus providing evidence 
against the formation of a Bragg glass in two dimensions. 

Further exploration of these disordered 2d lattice models with SOS rep- 
resentations will certainly help to address the fundamental issue of how 
non-random critical systems are affected by quenched disorder since most 
of these non-random models are critical[ll]. For example, the non-random 
FPL model flows to the densely-packed loop (DPL) fixed point of the O(n) 
model upon the perturbation of holes [ 17 1. How this DPL fixed point get 
modified by bond randomness (if at all) can now be examined by using non- 
bipartite matching algorithms. Another exciting extension of these polyno- 
mial algorithms to compute the energetics of dislocation loops in 3d random 
elastic media is now also feasible. Compared with the past few years in 
which polynomial and exact algorithms have been productively explored, it 
is fair to say that the next few years will see a rapid closing-in on a class of 
even NP-hard disordered systems which allow polynomial and approximant 
algorithms with near-optimal solutions of guaranteed bounds. 

We thank J. Kondev, C.L. Henley and A. A. Middleton for useful discus- 
sions. Part of this work is done in collaboration with D.S. Fisher which is 
also gratefully acknowledged. 
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